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Abstract. The partial differential equation governing the formation of the radiation spectrum in the 
shocked plasma accreting onto an x-ray pulsar is solved to obtain the associated Green's function describing 
the scattering of monochromatic radiation injected into the plasma at a fixed altitude above the star. 
Collisions between the photons and the infalling electrons cause both the ordered and random components 
of the plasma energy to be transferred to the radiation, which escapes by diffusing through the walls of the 
accretion column. The analytical solution for the Green's function provides important physical insight into 
the formation of the observed spectra in x-ray pulsars. Interesting mathematical aspects of the problem 
include the establishment of a new summation formula involving the Laguerre polynomials, based on the 
calculation of the photon number density via integration of the Green's function. 



I. INTRODUCTION 

The analysis of solutions to partial differential equations satisfying physically prescribed boundary 
conditions often leads to new insights into the relationships between the various special functions of 
mathematical physics. In this article, we employ the methods of classical analysis to obtain the closed- 
form solution for the Green's function describing the scattering of radiation inside the plasma accreting onto 
the surface of a rotating neutron star in an x-ray pulsar. The method is based on the separation of the 
two-dimensional partial differential equation into two second-order ordinary differential equations depending 
on the energy and space variables. By applying suitable boundary conditions in the spatial dimension, we 
are able to determine the eigenvalues and show that the corresponding spatial eigenfunctions are Laguerre 
polynomials. The solutions to the associated energy equation are Whittaker functions. The global solution 
for the two-dimensional Green's function is obtained by forming an expansion based on the energy and 
space eigenfunctions. The availability of this solution also allows us to calculate the distribution of the 
radiation escaping from the accreting gas, which forms the spectrum observed at earth. Beyond the direct 
physical relevance of the Green's function, the method of solution also yields several additional results of 
mathematical interest, including a new summation identity for the Laguerre polynomials. 

We provide a brief overview of the physical problem before proceeding with the main derivation. The 
fundamental power source for the radiation produced in bright x-ray pulsars is the gravitational accretion 
(inflow) of fully ionized plasma that is channeled onto the poles of the rotating neutron star by the strong 
magnetic field. In these luminous sources, the pressure of the photons governs the dynamical structure of 
the accretion flow, and therefore the gas must pass through a radiation-dominated shock on its way to the 
stellar surface. The kinetic energy of the gas is carried away by the high-energy radiation that escapes from 
the column, which allows the material to come to rest on the stellar surface. 1 ' 2 The strong compression of the 
infalling gas drives its temperatures up to several million Kelvins. The gas therefore radiates x-rays, which 
appear to pulsate due to the star's rotation. The observed x-ray spectrum is often distinctly nonthermal, 
indicating that nonequilibrium processes are playing an important role in the radiative transfer occurring 
inside the accretion column. 3 
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The nonthermal shape of the spectrum is due to the transfer of energy from the gas to the photons via 
electron scattering. The plasma possesses both "ordered" kinetic energy associated with the inflow, and 
"random" kinetic energy associated with the thermal motion of the particles. These two types of energy 
are transferred to the photons via first- and second-order Fermi processes, respectively. Our primary goal in 
this article is to obtain the exact solution for the Green's function describing the effect of electron scattering 
on monoenergetic seed photons injected from a source located at an arbitrary altitude inside the accretion 
column. The Green's function contains a complete representation of the fundamental physics governing the 
propagation of the photons in the physical and energy spaces. Since the transport equation governing the 
radiation spectrum is linear, we can compute the solution associated with any seed photon distribution via 
convolution. In x-ray pulsars, the most important sources of seed photons are bremsstrahlung, cyclotron, and 
blackbody emission. 4 The Green's function provides a direct means for investigating how photons produced 
by the various source mechanisms contribute to the observed nonthermal spectra in x-ray pulsars. 



II. FUNDAMENTAL EQUATIONS 



The plasma flowing onto the neutron star is assumed to have a cylindrical geometry, maintained by the 
strong dipole magnetic field. We define the spatial coordinate z as the altitude measured from the stellar 
surface along the axis of the cylindrical accretion column. The gas flows onto the star with velocity v(z) < 0, 
which vanishes at the surface of the star. The Green's function, f G (zo, z, eo, e), is defined as the radiation 
distribution measured at location z and energy e resulting from the injection of N photons per second with 
energy eo from a monochromatic source located at altitude zq inside the column. In a steady-state situation, 
/ G satisfies the Kompaneets transport equation 4 ' 5 



^ dfa = dve dfc_ d_ f c df G \ _ 
dz dz 3 de dz \ 3n P a\\ dz J i rsr 



n e ac 1 d 
m e c 2 e 2 de 



eUf a +kTr Ja 
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where r, is the radius of the accretion column, i osc represents the mean time photons spend in the plasma 
before escaping through the walls of the column, cf denotes the angle-averaged electron scattering cross 
section, c is the speed of light, k is Boltzmann's constant, and T e , n e , and m e denote the electron 
temperature, number density, and mass, respectively. The left-hand side of ([I]) represents the comoving 
time derivative of the Green's function, and the terms on the right-hand side represent first-order Fermi 
energization ("bulk Comptonization" ) , spatial diffusion along the column axis, photon escape, stochastic 
(thermal) Comptonization, and photon injection, respectively. 

The transport equation employed here is similar to the one analyzed by Becker 1 and Becker and Wolff, 2 
except that thermal Comptonization has now been included via the appearance of the Kompaneets (1957) 
operator, which is the fourth term on the right-hand side of ([1}. The radiation number density n G and 
energy density U G associated with the Green's function are given by 

/>oo />oo 

n G (z)= e 2 f G (z ,z,e ,e)de , U G {z) = I e 3 f G (z , z, e , e) de . (2) 
Jo Jo 

The Green's function provides fundamental physical insight into the spectral redistribution process, and it 
also allows us to calculate the particular solution for the spectrum f(z, e) associated with an arbitrary photon 
source Q(z, e) using the integral convolution 6 

f(z, e) = r r /o(wo,e) el Q(z ,eo) de, dz, , (3) 
Jo Jo iv 

where the source function Q is normalized so that e 2 Q(z, e) dedz gives the number of seed photons injected 
per unit time in the altitude range between z and z + dz and the energy range between e and e + de. 

Following Lyubarskii and Sunyaev, 7 we will assume that the electron temperature T e has a constant value, 
which is physically reasonable since most of the Comptonization occurs in a relatively compact region near 
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the base of the accretion column. In this situation, it is convenient to work in terms of the dimensionlcss 
energy variable Xi defined by 

X(e) = , (4) 
and the dimensionless optical depth variable t, defined by 

dr = n e (z) crii dz , r(z) = / n e (z') an dz' , (5) 

Jo 

where z and r both vanish at the stellar surface. The electron number density n e is related to the inflow 
velocity v via 

M = tttq m p n e \v\ = constant , (6) 

where M is the mass accretion rate onto the magnetic pole of the star and m p is the proton mass. Note 
that we can write the Green's function as either f G (zq, z, cq, e) or / g (tq, t, XQj x) since the variables (z,e) 
and (Xj t ) are interchangeable by utilizing (TjE]) and (JSJ). 

Making the change of variable from z to r in |T]), we find after some algebra that the transport equation 
for the Green's function can be written in the form 



vdfe = l_dvxdfg + 1 d 2 f G 

c dr c dr 3 dx 3 dr 2 c 2 



fa 



a kT e 1 d 



N S{x- Xo)S(t ~ t )_ ^ 



irr 2 ckT e e 2 



au m e c 2 x 2 dx 

where to = t(zo), xo = x( e o), and we have introduced the dimensionless constant 

_ irr m p c 

* = M{ n <, ± )V* ' (8) 

which determines the importance of the escape of photons from the accretion column. It can be shown that £ 
is roughly equal to the ratio of the dynamical (accretion) timescale divided by the timescale for the photons 
to diffuse through the column walls, 4 and therefore the condition £ ~ 1 must be satisfied in order to ensure 
that the radiation pressure decelerates the gas to rest at the stellar surface, with the kinetic energy of the 
material carried away by the photons escaping through the column walls. 8 

Following Becker and Wolff, 2 we compute the mean escape time using the diffusive prescription 

icsc(z) = ^— - , t±(z) = n e a ± r , (9) 
c 

where tj_ represents the perpendicular scattering optical thickness of the cylindrical accretion column. Note 
that tj_ and t esc are each functions of z through their dependence on the electron number density n e (z). 
Becker 8 confirmed that the diffusion approximation employed in ([9]) is valid since t± > 1 for typical x-ray 
pulsar parameters. The technical approach used to solve for the Green's function, carried out in § 3, involves 
the derivation of eigenvalues and associated eigenfunctions based on the set of spatial boundary conditions 
for the problem. 9 ' 10,11 ' 12 



III. EXACT SOLUTION FOR THE GREEN'S FUNCTION 



Lyubarskii and Sunyaev demonstrated that the transport equation Q is separable in energy and space 
if X 7^ Xo an d the velocity profile has the particular form 

v(t) = —acT , (10) 

where a is a positive constant. This expression provides a reasonable approximation of the actual velocity 
profile in an x-ray pulsar accretion column. By combining ([5]), ([S]), ([5]), and (jlOp . we can express r as an 
explicit function of the altitude z, obtaining 



1/4 / o . \ 1/2 
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Using this result to substitute for r in (flT))) . we note that the velocity profile required for separability is 
related to z via 



v(z) = - ~ 



1/4 /0 \ 1/2 



(12) 



This profile describes a flow that stagnates at the stellar surface (r = 0, z = 0) as required. 

Adopting the velocity profile given by (fTU|) , we find that the transport equation ([7]) for the Green's function 
can be reorganized to obtain 



a X 9f G a kT e 1 d 
3 dx cri m e c 2 x 2 9x 



1 d 2 f df 
3 9r 2 9r 



■^Vo5(x-Xo)^(t-to) 



7rr 2 cfcT e eg 



(13) 



When x Xo, the (^-function source term on the right-hand side makes no contribution, and the differential 
equation is therefore linear and homogeneous. In this case the transport equation can be separated in energy 
and space using the functions 

f\(T,x)=g(\T)h(\x) , (u) 

where A is the separation constant. We find that the spatial and energy functions, g and h, respectively, 
satisfy the differential equations 

" " ' ' (15) 



1 d 2 g dg 



- e^r 1 )g = , 



and 



_L__! 

X 2 dx 

where the parameter [3 is defined by 



dx 



dx 



(16) 



., - 1T - (17) 

It can be shown that (3 determines the relative importance of bulk and thermal Comptonization. When is 
of order unity, the two processes are comparable, and when j3 > 1, the bulk process dominates. 



3.1 Eigenvalues and Spatial Eigenfunctions 



In order to obtain the global solution for the spatial separation function g, we need to understand the 
physical boundary conditions that g must satisfy. In the downstream region, as the gas approaches the stellar 
surface, we require that the advective and diffusive components of the radiation flux must both vanish due to 
the divergence of the electron density as the gas settles onto the star. The advective flux is indeed negligible 
at the stellar surface since v — ► as r — > [see Eq. (flTj|) ]. However, in order to ensure that the diffusive 
flux vanishes, we must require that dg/dr — > as r — * 0. Conversely, in the upstream region, we expect 
that g — > as r — > oo since no photons can diffuse to large distances in the direction opposing the plasma 
flow. With these boundary conditions taken into consideration, we find that the fundamental solution for 
the spatial separation function g has the general form 



i(X,t) oc 



' e - a (3+ W )r 2 /4 V ^ 1 ) ; T > ^ 

e ^(3+-)T 2 /4 M ( p; 1,2-d) > T < ro 



where M and U denote confluent hypergeometric functions, 13 and we have made the definitions 



2A 
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(9 + 12^ 



,1/2 



(18) 



(19) 



Equation (|15[) is linear, second-order, and homogeneous, and consequently both the function g and its 
derivative dg/dr must be continuous at the source location, t = tq. The smooth merger of the M and U 
functions at the source location requires that their Wronskian, uj(t), must vanish at r = to, where 



1 awr 



1 awr 2 



1 awr 



1 awr 
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This condition can be used to solve for the eigenvalues of the separation constant A. By employing 
equation (13.1.22) from Abramowitz and Stegun 13 to evaluate the Wronskian, we obtain the eigenvalue 
equation 

r(|)(2^ eWT2/2 = 0j 
r (p) 

The left-hand side vanishes when T(p) — > ±00, which implies that p = —n, where n = 0, 1, 2,. . . By combining 
this result with (|19p . we conclude that the eigenvalues A„ are given by 

Anw + w + 3 . . 

A„ = , n = 0,l,2,... (22) 

When A = A„ , the spatial separation functions g reduce to a set of global eigenfunctions g n that satisfy the 
boundary conditions at large and small values of r. In this case, we can use equations (13.6.9) and (13.6.27) 
from Abramowitz and Stegun 13 to show that the confluent hypergeometric functions M and U appearing in 
(|20p are proportional to the generalized Laguerre polynomials L„ X ^ 2 \ and consequently the global solutions 

awr 2 



for the spatial eigenfunctions can be written as 

(23) 



3 „(r)^ 5 (A„,r)=e^ 3 +^ 2 / 4 4- 1/2) ( 



Based on equation (7.414.3) from Gradshteyn and Ryzhik, 14 we note that the spatial eigenfunctions satisfy 
the orthogonality relation 



r(n+l/2) 

— — n = m 



e 3ar/2 9n (T)g m (T)dT = {f^> (24) 

10, n f= m . 



3.2 Energy Eigenfunctions 



The solution for the energy separation function h depends on the boundary conditions imposed in the 
energy space. As x ~* 0, we require that h not increase faster than x -3 since the Green's function must 
possess a finite total photon number density [see Eq. l[2"])]. Conversely, as x ~ * 00 1 we require that h decrease 
more rapidly than x~ A m order to ensure that the Green's function contains a finite total photon energy 
density. Furthermore, in order to avoid an infinite diffusive flux in the energy space at x — Xo> the function 
h must be continuous there. The fundamental solution for the energy eigenfunction that satisfies the various 
boundary and continuity conditions can be written as 

h n {x) = h{X n ,x) = < K _ 4 _ W2 ,, , , ^ ( 25 ) 

where M„ ijU and W K<I1 denote Whittaker's functions, and we have made the definitions 

K =I(/? + 4), [(3-/3) 2 +4/3A„] 1/2 . (26) 

Note that each of the eigenvalues A„ results in a different value for /x, and the parameter (3 is defined in (| 1 T|) . 
Equation (I25p can also be written in the more compact form 

K{X) = X K ~ 4 e- x/2 M K , M ( Xm m) W KtlA ( Xmax ) , (27) 

where 

Xmin = min(x, Xo) , Xmax = max(x, Xo) ■ (28) 
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3.3 Eigenfunction Expansion 

The spatial eigenfunctions g n (y) form an orthogonal set, as expected since this is a standard Sturm- 
Liouville problem. The solution for the Green's function can therefore be expressed as the infinite series 

oo 

X ) = ]T C n9n {T)h n { X ) , (29) 

n=0 

where the expansion coefficients C n are computed by employing the orthogonality of the eigenfunctions, 
along with the derivative jump condition 



lim -fS. 



dfc 
dx 



3N pkT e S(T-T ) 

2 4 ' l dU J 



a tttq c Cq 



X=Xo+e 

which is obtained by integrating the transport equation () X3[) with respect to \ m a small region surrounding 
the injection energy \o- Substituting using the expansion for f G yields 



lim r, Yl Cn 9n(r) [h' n (xo + e) - h'Jxo - s)] = - '» ^ "-- y v ^ — — , (31) 



3N PkT e 6(T-T ) 



n=0 



where primes denote differentiation with respect to x- By employing (|27j) for h n , we find that 



n=0 



(32) 



where the Wronskian, 20, is defined by 

an(xo) ee m k m w^( X o) - w K ,M m' k j X o) ■ (33) 

We can evaluate 2U(xo) analytically using equation (2.4.27) from Slater 15 , which yields 

r(i + 2p) 

2»(Xo) = - =7 —TV ■ (34) 

Using this result to substitute for 2U(xo) m H32|) and reorganizing the terms, we obtain 

^ r(f + 2^) C n9n (T) 3N (3e^/ 2 S(T-T ) 

2_, I>- K+ i) aTr^cXoMW ' > 

where fi is a function of A n via ((26]). We can now calculate the expansion coefficients C n by utilizing 
the orthogonality of the spatial eigenfunctions g n represented by (|24|) . Multiplying both sides of (|35|) by 
e 3ar / 2 (? m (T) and integrating with respect to r from zero to infinity yields, after some algebra, 

= 3iV /3V2^e*°/2 e ^/2 r ( M _ K + i) n! g „( To ) 

7rr 2 c(fcT e )3xSV^ r(l + 2/x) r(n + |) ' ( ' 

The final closed- form solution for the Green's function, obtained by combining (|2"?| . (12T))) . and (|3"6"|) , is given 

by 

3iV /3e 3QT o/ 2 V2^X K ^ 4 e (xo ~ x)/2 r(/i-/t+±)n! 



/ G (To,r,xo,x) - nr*c(kT e f x ^ ^ r(l + 2 M ) T(n + ±) 

X 5n(T-o)5n(T)M KiA1 (Xmin) W K , M (Xmax) , (37) 

where the spatial eigenfunctions g n are computed using (|23[) and the parameters k and are given by (|26[) . 
This exact, analytical solution for / G provides a very efficient means for computing the steady-state Green's 
function resulting from the continual injection of monochromatic seed photons from a source at an arbitrary 
location inside the accretion column. The eigenfunction expansion converges rapidly and therefore we can 
generally obtain an accuracy of at least four significant figures in our calculations of / G by terminating the 
series after the first 5-fO terms. 
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Figure 1. Green's function X /g describing the photon spectrum inside an x-ray pulsar 
accretion column [Eq. (|37[) ] plotted as a function of the dimensionless photon energy x f° r 
iV = 1, T e = 10 7 K, r = 10 4 cm, a = 0.1, £ = 1.5, xo = 0.1, r = 0.5, and (a) = 0.3, 
(b) (3 — 1.5. The values of the optical depth are r = 0.01 (solid line), r = 1.0 (dashed line), 
and r — 1.5 (dot-dashed line). 




Figure 2. Same as Fig. 1, except xo = 0.5. In this case the energy of the injected photons 
is comparable to the thermal energy of the electrons. 



3.5 Numerical Examples 



In this section we illustrate the computational method by examining the dependence of the Green's 
function f G (to, t, xo> x) 011 the optical depth r and the dimensionless energy \. We remind the reader that 
the solution for the Green's function represents the photon spectrum inside the accretion column at the 
specified position and energy, resulting from the injection of monochromatic photons with dimensionless 
energy xo from a source located at optical depth To. Analysis of / G therefore allows us to explore the effects 
of bulk and thermal Comptonization as the photons diffuse throughout the plasma, eventually escaping 
through the walls of the accretion column. In terms of x an d t, the radiation number density associated 
with the Green's function is given by [cf. Eq. 

/>oo 

n G (r) = (kT e f / x 2 f G (ro,T,xo,x)d X ■ (38) 
Jo 

It follows that (kT e ) 3 x 2 f G d-X equals the number of photons per unit volume at optical depth t with 
dimensionless energy between x and X + ^X- 

In Fig. 1 we use (|37|) to plot x 2 / G as a function of x an d r for the parameter values No = 1, T e = 10 7 K, 
ro = 10 4 cm, a = 0.1, £ = 1.5, xo — 0.1, to = 0.5, and either j3 = 0.3 or (3 — 1.5. We remind the reader 
that the value of (3 determines the relative importance of bulk and thermal Comptonization [see Eq. (|17|) ]. 
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Both spectra in Fig. 1 display a peak at the energy of the injected photons, x — 0-1. F° r the case with 
(3 — 0.3, there is also a noticeable Wien hump in the spectrum around x ~ 2 due to the effect of thermal 
Comptonization. 3 However, when {3 = 1.5, the Wien feature is completely hidden by the power-law shape 
associated with the dominant bulk Comptonization process. The Green's function is plotted in Fig. 2 for the 
same parameters used in Fig. 1, except that we now set xo = 0.5, so that the energy of the injected photons, 
eq, is comparable to the thermal energy of the electrons, kT e . In this case, the Wien hump is still visible 
when (3 — 0.3, although it has begun to merge into the peak associated with the injected photons because 
of the larger value of xo- Note that all of the spectra in Figs. 1 and 2 display a normalization that decreases 
with increasing r, due to the fact that the downward advection of the electrons tends to "trap" the photons 
near the bottom of the accretion column, close to the stellar surface. 8 

Equation (f3"T|) can also be used to calculate the spectrum of the radiation escaping through the walls of the 
accretion column, which comprises the radiation spectrum observed at Earth. Because the x-ray pulsars are 
distant cosmic sources, the observed spectra are the result of emission escaping over the entire vertical extent 
of the accretion column, and therefore we must perform a vertical integration with respect to r in order to 
compare the model spectra with actual x-ray observations. This procedure is carried out in reference [4], 
where the resulting spectra are compared with the x-ray data for several astrophysical sources, confirming 
excellent agreement. 



IV. IDENTITIES INVOLVING THE LAGUERRE POLYNOMIALS 



We are primarily interested here in the mathematical properties of the solution for the Green's function 
given by (|37p . In the present application, it is possible to compute the photon number density, n G , either by 
integrating the series expansion for the Green's function term-by-term, or by solving directly the ordinary 
differential equation satisfied by n G . By equating the results obtained for the photon number density using the 
two methods, we can derive an interesting new summation identity for the Laguerre polynomials ll l ^ (x) 
appearing in the expression for the spatial eigenfunctions g n (T~). We carry out this procedure below. 



4.1 Integration of the Green's Function 



By using (|37p to substitute for f G in (|38[) and integrating term-by-term, we find that the associated photon 
number density n G can be written as 

"° (t) = £ r(l + 2M)r(n + l) M 9n {T)K M , (39) 

where 

rXo roo 
K G ( X o) = W K Mo) x K - 2 e-* /2 M K M)dx + M K Mo) X^ 2 e~ x/2 W Kjfi ( X ) d X • (40) 



The two indefinite integrals on the right-hand side of (j4H)) can be evaluated using equations (3.2.6) and 
(3.2.12) from Slater 15 . After some algebra, we find that 

Xo p-Xo/2 v «-l 

X K - 2 e- x/2 M K M) dX = , X \ M K .^( X o) , (41) 

K + fl — 5 

and 

X K ~ 2 e~ x/2 W K M) dX = e-* o/2 Xo' 1 W-i^Oeo) . (42) 

Xo 

Combining relations yields 

e -Xo/2 y K-l 

^g(xo) = , % W K ,MM K -iAxo) + er Xo/2 xt M K Mo)W K -x,M , (43) 

K + /J, — o 
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or, equivalently, 

^g(xo) = ^- r [K-iMo) W K -iAxo) - <-i lM (xo) M„_ 1)/t (x ) , (44) 

where we have used equations (13.4.32) and (13.4.33) from Abramowitz and Stcgun, 13 and primes denote 
differentiation with respect to \o ■ The Wronskian on the right-hand side of (|44|) can be evaluated using (|34|) 
to obtain 

^^^Hr^^D' (45) 

which can be further simplified by applying (j2Ul) . yielding 

xge-*°/2 r(i + 2 M ) 

^g(Xo) = -T7T oT ^7 — TT . (46) 

where the eigenvalues A„ are given by (12"2")1 . By using (|4"b]l to substitute for K G (xo) in (|3"5|) and simplifying, 
we find that the photon number density can be evaluated using the expansion 

_ 3N e 3 ^ 2 o/ 2 V2^ ^ nlg n (T )g n (r) 
Hg[T) KTgc^a ^ r(n+|)(A„-3) ■ 10 

This expression allows the computation of the photon number density at any optical depth t inside the 
accretion column for given values of the photon injection rate Nq and the source location tq. 



4.2 Solution of the Differential Equation 

The availability of the transport equation (|13[) provides us with an alternative means for computing the 
photon number density, n G , by directly solving the differential equation 

1 r,2n dn r , ^999 NqS(t~t ) , . 

aT—^ + an G -£ 2 a 2 T 2 n G = — ^ ^ , (48) 



3 dr 2 dr G * G irr 2 c 

which is derived by operating on (flU|) with (kT e ) 3 J °° x 2 dx- The homogeneous version of (|4"8"|) obtained 
when t To admits the fundamental solutions 

|4e-^%(a,|,^) , r>r , 

where A and B are constants and 

w — 3 _ , 2 n 1/2 



a — 



- S (9+120 V . (50) 

The solutions in (|4"5)) are consistent with the physical boundary conditions at large and small r discussed in 
§ 3.1 [sec Eq. (p~5|) ] . The constants A and -B are determined by requiring that n G (r) be continuous at the 
source location tq, and that it satisfy the derivative jump condition 



dn r , 
hm — - 

e^o dr 



dn G 
~ ~~dT 

T=T +e 



^ , (51) 

TTTq c 



obtained by integrating (j4"8")) with respect to t in a small region around tq. Combining (|4U|) and (j5Tj) with 
the continuity condition and the expression for the Wronskian given by (|2ip , we find after some algebra that 



A = f°^ e^oW a ,l, ^f] , (52) 

7T7QCV 2,TtOtW \ £ £ J 

B = f°^ e ^o/ 4 L/fa,l,^i) , (53) 
7rrQCv 27raiiJ V z z / 
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and consequently the global solution for the photon number density is given by 

= ZNoTja) eQ(3 -. )ToV 4 e - Q (3 + ^/4 M (12^) u L^ZEL*) (54) 
irrficV2naw V. 2 2 / V 2 2 / 

where 

r min = min(r, r ) , r max = max(r, t ) . (55) 



4.3 Summation Identity for the Laguerre Polynomials 



Equations (|47|) and (f54|) provide two independent expressions that can each be used to calculate the 
photon number density n G . This fact allows us to derive an interesting new summation identity involving 
the Laguerre polynomials L„ 1 ^ 2 '(x) appearing in equation (|23p for the spatial eigenfunctions g n . By setting 
(|47f and (|54|) equal and substituting for A„ using (f22|) . we obtain after some simplification 



(-1/2) 



(ssp) _ r(a) 



n=0 



(n + o)T(n+ |) 



1 aror: 



Setting x = awr 2 /2 and = awTg/2 yields the more compact form 



Jll l^n 

(^+"«)r(n + i) 



(a;o)4~ 1/2) (a0 _ r(a) 



M | a, -, x r 



(56) 



where 

a: m in = min(iE, x ) , x max = max(a;, x ) . 

Equation ([57)) has not appeared in the previous literature, although Exton, 17 Srivastava, 18 and Manocha 19 
have derived several related identities. 



(57) 
(58) 



V. CONCLUSION 

In this article we have applied the methods of classical analysis to derive the closed-form solution for 
the Green's function, f G , given by (|3"T|) , describing the bulk and thermal Comptonization of monochromatic 
seed photons scattered by hot, infalling electrons inside an x-ray pulsar accretion column. X-ray pulsars are 
rotating neutron stars, which are the most dense and compact solid objects known to exist in the universe. 
These enigmatic sources are characterized by super strong magnetic, gravitational, and radiation fields, 
making them the most extreme physical "laboratories" in the universe. It is therefore of great theoretical 
interest to obtain the best possible understanding of the spectral formation process in x-ray pulsars, because 
the observed radiation provides us with the only available window into their physical nature. As demonstrated 
in Figs. 1 and 2, the Green's function is characterized by a power-law shape at moderate photon energies, 
with an exponential turnover at higher energies, in agreement with the spectra of many sources. 4 We conclude 
that bulk and thermal Comptonization in the shocked pulsar accretion column provides a natural explanation 
for the typical high-energy spectra produced by x-ray pulsars. 

Our primary interest here is in exploring the mathematical properties of the Green's function. In 
particular, we established that by performing two independent calculations of the photon number density 
n G associated with / G , one can obtain an interesting new summation identity involving the Laguerre 

polynomials, L„ 1 ^ 2 \x), expressed by (|57|) . The derivation is based on the simultaneous calculation of 
n G using either term-by-term integration of the Green's function expansion (|37p . which yields (|47[) . or 
instead via the direct solution of the differential equation for n G given by AS)) , which leads to (|54[) . The 
new identity for the Laguerre polynomials obtained here is related to various similar expressions presented 
by Chatterjea, 16 Exton, 17 Srivastava, 18 Manocha, 19 Varma, 20 and Gradshteyn and Ryzhik, 14 although our 
results are distinct from theirs. We expect that our new expression may be of potential benefit in various 
problems of mathematical physics. 

The author would like to gratefully acknowledge the helpful comments provided by the anonymous referee. 
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